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ABSTRACT 

This  technical  report  analyzes  the  performance  of  various 
methods  of  Squared  Magnitude  Coherency  (SMC)  estimation  when  the 
available  data  is  very  short. 

Three  methods  of  Squared  Magnitude  Coherency  estimation  were  exam¬ 
ined.  The  first,  Coherency  estimation  using  the  periodogram  and  a  FFT 
(the  frequency  domain).  The  second,  Fitting  an  AR  model  to  the  process, 
estimanng  its  coefficients,  and  calculating  the  SMC  from  the  coefficients 
(the  time  domain  analysis).  The  third  method  is  a  direct  coherency  esti¬ 
mation  via  constrained  least  squares  linear  prediction  proposed  by  A.  H. 
Nuttall. 

For  a  short  data  the  SMC  estimation  using  the  time  domain  method  has  a 
large  bias  error,  and  a  large  confidence  interval.  The  other  two  methods 
have  reasonable  performance  for  short  data. 

Most  of  the  previous  work  that  calculated  the  SMC  estimation  confidence 
interval,  assumed  that  to  the  estimated  SMC  is  Gaussian  distributed  (  Ben- 
dat  &  Piersol,  Brockwell  &  Davis  ).  Simulation  results  show  that  this 
assumption  holds  only  for  SMC  values  between  0.3  -  0.8,  for  long  data  as 
well  as  for  short  data. 

Simulation  results  reveal  that  as  soon  as  the  number  of  samples  is  such 
that  the  law  of  large  numbers  holds  (N  greater  then  =  40),  all  the  expres¬ 
sions  of  the  estimated  SMC  density  function  developed  by  Carter  et  al. 
and  Sakai  et  al.  still  hold. 
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INTRODUCTION 


Squared  Magnitude  Coherence  (SMC)  and  time  delay  are  two  quantities  needed 
for  passive  detection  and  localization  of  a  radiating  source.  The  source  signal,  together 
with  environmental  noise,  is  received  by  sensors  at  known  locations.  The  coherence  func¬ 
tion  is  useful  for  detection  and  sepiU'ation  of  signals,  since  it  measures  the  causality  and 
similai'ity  between  a  pair  of  signals  received  by  two  different  sensors. 

In  the  case  of  two  channels,  a  ver^'  common  test  is  to  compare  the  SMC  estimate 
computed  using  sample  sequences  from  the  two  channels  to  a  threshold^.  This  is  an  espe¬ 
cially  useful  test  because  the  statistical  behavior  of  the  SMC  estimate  is  well  understood 
in  the  absence  of  a  common  signal,  allowing  thresholds  coiresponding  to  particular  false 
alarm  probabilities  to  be  readily  established. 

When  designing  an  underwater  acoustic  telemetry  system  using  a  chirp  D^SK 
modulation  scheme,  the  choherence  between  two  received  pulses  is  very  important  in 
order  to  evaluate  the  system  performance  and  gives  a  great  deal  of  information  on  the 
channel  characteristics.  As  only  a  short  data  interval  is  available  (each  bit  is  veiy  short),  a 
question  arising  regarding  the  quality  (in  the  sense  of  bias  and  confidence)  of  the  conven¬ 
tional  SMC  estimation  methods,  and  whether  then  performance  analysis  for  short  data 
intervals  still  holds. 

Given  two  stationary  zero  mean  random  processes  .v(0  and  y(/),  the  squared 
magnitude  coherence  (SMC)  is  defined  as  : 


lC((o)r-= 


(1.1 


Where  .V^to))  is  the  cross  power  spectrum  between  the  signals  .v{/)  and  yto,  and  S'crlw) 
and  Syy(w)  are  the  auto  power  spectra  of  .v(/)  and  y(i)  respectively.  The  most  traditional 
way  of  SMC  estimation  uses  the  peiiodogram  and  compute  the  individual  spectra  usually 
via  the  Fast  Fourier  Transform 


I 


All  of  the  methods  based  on  estimation  of  the  auto  and  cross-spectra  have  'he 
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advantage  of  not  assuming  any  signal  model,  i.e  no  a  priori  information  about  the  signal  is 
needed.  The  disadvantage  of  these  methods  is  that  in  cases  where  there  only  is  limited 
data  length,  a  tight  tradeoff  between  resolution  and  estimate  stability  exists. 

A  new  modeling  approach  to  spectral  estimation  was  proposed  in  the  late  sixties 
by  Parzen  and  Burg^^.  This  new  method  has  based  on  the  idea  that  there  is  always  any  a 
priori  knowledge  about  the  processed  signals,  and  uses  this  information  to  fit  a  mathemati¬ 
cal  model  to  the  process.  The  drawback  of  this  approach  is  the  need  to  fit  the  best  model 
and  the  needto  find  the  closest  model  order. 

Many  deterministic  and  stochastic  discrete  time  processes  encountered  in  practice 
are  well  approximated  by  a  rational  transfer  function  model.  In  this  model  an  input  driving 
sequence  {  n„  >  and  the  coiTesponding  output  sequence  {  .v„  >  that  is  to  model  the  data 
ai’e  related  by  the  linear  difference  equation, 

^  Ibi  1/r,.,-  - ^  [aj ( 1 .2) 
where  x  and  n  are  [d  x  l]  vectors,  and  the  [a]’s  and  [b]’s  are  [d  .r  d]  matrices. 

The  cross  specti'al  mamix  follows  from  Eq.  (1.2); 

[S(f)]=[A(f)r'iD(f)r  iz][B{f)v  iA(f)r'  (1.3) 

where: 

[A  (/  )I  =  [/ 1  +  ^  [o  It  exp  (  -jl-nkf  ) 

[B(f)]  =  U]  +  '^  lb]kC\p(-j2Kkf  ) 
and 

IE]  is  the  driving  sequence  covariance  matrix. 

Methods  of  fitting  a  model  and  its  order  using  measured  data  exist  Kay^^  (page 
1395)  Wax^^  and  Marple^^.  In  the  most  common  models  it  is  assumed  that  the  received 
time  series  can  be  described  as  an  autoregressive  (AR)  process,  moving  average  (MA)  pro¬ 
cess,  or  combination  of  them  (autoregressive  moving  average  (ARMA)  procesx)^’ 
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Under  the  assumption  of  fitting  a  linear  model  to  the  measured  process,  one 
method  of  estimating  the  SMC  is  by  estimating  the  cross  spectral  matrix  and  substituting 
the  estimated  values  in  Eq.  (1.1  Some  "semi  direct"  methods  of  estimating  the  SMC 
based  on  AR,  MA  or  ARMA  process  modeling  were  suggested  by  Dae  Hee  Youn  at  el.^"^, 

and  Nuttall^^.  In  these  papers,  methods  of  directly  estimating  the  ratio  and  \ 

were  suggested. 

The  estimator  statistical  properties  ai'e  give.i  for  only  a  veiy  few  estimation 
methods.  The  most  detailed  statistical  analysis  is  given  for  the  SMC  estimation  technique 
using  the  FFT  for  the  cross  spectral  matiix  estimation.  The  probability  function  of  the 
SMC  estimator  and  a  numerical  method  how  to  calculate  it  are  described  in  papers  by 
Carter  at  el^'^*^^.  The  bias  of  the  estimator  is  calculated  in  and  the  confidence  inter¬ 
val  is  analyzed  by  Bendat  &  Piersol^  (chap  6),  by  Fay^  by  Carter^^  and  by  Brockwell  & 
Davis^.  The  statistical  properties  of  Squtmed  Magnitude  Coherence  estimation  via  Mul¬ 
tivariate  AR  modeling  is  given  by  Sakai  &  Tokumaru^®. 

Two  important  assumption  were  made  in  all  the  statistical  analyses  of  the  SMC 
estimator.  The  first  assumption  is  that  the  available  data  is  very  long  such  that  the  central 
limit  theorem  holds,  and  the  second  being  that  the  driving  signal  is  white  Gaussian  noise 
(WGN) 

In  many  applications  such  as  under  w'ater  explosion,  seismic  signals,  or  other  man- 
sient  events,  the  received  data  for  analysis  is  very  short  by  definition  and  the  confidence 
interv'al  using  the  traditional  FFT  method  is  veiy  huge.  In  this  report,  we  summarize  the 
work  done  in  examining  the  "performance"  of  methods  of  SMC  estimation  using  modem 
spectral  estimation  techniques  and  find  their  limitation  when  only  short  data  is  available. 
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H  STATISTICAL  METHODS  OF  SMC  ESTIMATION 

II.l  SMC  ESTIMATION  VIA  THE  PERIODOGRAM 

E.  J.  Hannan  proved  that  for  a  time  series  x(n)  which  is  generated  by  a  linear 
process  and  has  a  continuous  spectrum  which  is  not  zero  at  /  ,  the  asymptotic  (the  total 
number  of  data  points  is  going  to  infinity)  disuibution  of  the  SMC  estimator  has  the  den¬ 
sity: 


(2.1) 

for  f  ^Q,n 

:  P  =  1, 

m  =  -j 

for  /  =0, 

:  P  =  y. 

-  = 

for  f  =n 

■  P  =  i' 

m=^N 

Here,  is  the  beta  function  and  iFi  is  the  confluent  hypergeometric  function,  N  is  the 
total  number  of  data  points. 

When  using  the  FFT  technique  as  a  SMC  estimator,  each  time  series  is  paititioned 
into  nj  equal  length  segments,  each  having  n  data  points.  The  segments  may  be  overlap¬ 
ping  or  disjoint  and  each  segment  is  multiplied  by  a  window.  The  FFT  of  the  weighted  « 
points  is  used  to  estimate  the  auto  and  cross  power  spectral  densities  which  ai'e  used  to 
form  the  SMC  estimate. 

for  rtd  defined  as  :  tij  =  ~  ,  the  density  function  of  the  SMC  estimator  is  ^  : 

-  1)  (1 -cYmI  -  ;  0=<cV^<l  (2.2) 

It  is  important  to  emphasis  that  Eq.  (2.1)  and  (2.2)  hold  only  in  the  limit  when  the  number 
of  samples  goes  to  infinity^  For  all  practical  purposes  these  expressions  can  be  used 
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when  the  number  of  points  is  large  enough  such  that  the  law  of  large  numbers  holds. 

Bendat  &  Piersol  ^  give  an  expression  for  the  confidence  inteiwal  when  the  FFT  method  is 
used  for  estimating  the  SMC.  In  their  analysis  they  assumed  that  the  estimated  SMC  is 
noiTnally  distributed^  [  page  312  in  the  second  ed  ],  an  assumption  that  holds  only  for 
values  of  SMC  approximately  between  0.35  and  0.95  and  for  mj>40^.  [  page  193  in  the  Is^ 
ed.  ]. 

Scannel  &  Cai'ter^^  gave  a  more  accurate  confidence  inteiwal  to  the  SMC  estima¬ 
tor  as  a  function  of  nu  and  the  confidence  level  a  .  If  the  upper  and  the  lower  confidence 
interval  are  defined  as  Al(C)  and  y4„(C)  respectively,  then 

probiAiiC)  <C  <A„(C)|C]  =  a 

where  a  is  the  desired  confidence  level.  They  found  this  inteiwal  numerically  by  taking  the 
estimated  SMC  density  function  as  in  Eq,  (2.2).  The  following  figures  give  their  results 
for  confidence  level  of  80%  and  95%  and  forn^  =  8,16,32,64  and  128. 
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II.2  STATISTICAL  PROPERTIES  OF  MULTIVARIATE 
AUTOREGRESIVE  SPECTRAL  ANALYSIS 


A  great  deal  of  effort  has  been  directed  toward  finding  a  SMC  estimator  using 
modern  power  estimation  method  having  better  performance  (in  the  sense  of  confidence 
inteiwal  and  frequency  resolution)  for  a  given  length  of  data.  This  problem  arises  espec- 
islly  in  cases  when  only  short  data  is  .available. 

One  way  of  solving  this  problem  is  to  model  the  analyzed  process.  The  most  gen¬ 
eral  linear  model  is  the  autoregresive  moving  average  (ARMA)  model. 

Using  an  ARMA  model  to  describe  the  process  is  problematic  since  a  set  of  nonlinear 
equiations  is  required  to  be  colved  in  order  to  estimate  the  coefficients  of  the  model  from 
the  mea>  jred  data.  It  is  easy  to  show  that  any  ARMA  process  can  be  modeled  as  an 
autoregresive  (AR)  process  with  an  infinite  number  of  coefficients.  Taking  a  sufficient 
number  of  coefficients  enable  us  to  use  an  AR  model  for  an  ARMA  process. 

Let  us  consider  a  two  dimensional  p-th  order  Gaussian  AR  process: 

|^[ayl.v„_^,  =  I/rl/f„  ^2.3) 

where  [ool  and  [i>l  are  the  unit  matices  ([/]),  and  [O;]  are  I?-r2]  AR  coefficients  matrices, 
and  {  }  is  a  sequence  of  Gaussian  white  noise  [Ivl]  vectors  with  mean  0  and  covari¬ 

ance  matrix  [I]. 

The  estimated  coefficient  matrices  satisfy  the  following  Yule  -  Walker  equation: 


(ail.  [i^i] . la„] 


/to  /to 
/t.,  /t-, 


/t  -p  +  1  /t  -^  + 1 


/to 

/t-i 


li-p.x 


RuRi . Rp 


(2.4) 


for  simplicity,  let  us  write  Etj.  (2.4)  as: 
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[ailJtfz] . ,[d,,]  [f]  =  y 


where 


Ro  Rc 
R-,  R- 


[n=  . 


R-o*i  R 


•p+\ 


Rj  is  defined  as  • 


y=  RuRi, . ,Rp 


1  1  _  _  j- 


The  estimated  spectral  density  matrix  [5V)]  of  Eq.  (2.3)  is  given  by  : 


[S{f)]  =  ([A (/)]')-'  [I]  (IAV)]-’)^ 


where; 


[A(f)]  =  ll]  +  ^  [d]t  exp {-jlKf  ) 


lt]  =  R^  +  ^ 


At  R-t 


The  SMC  estimator  is; 


izif  ) 

d\\{f)f  22(f  ) 


where;  %  is  the  ij  element  of  [S(f)] 


Following  Sakai  &  Tokumaru^^-^*,  the  SMC  estimator  distribution 


given  by  ; 


Fid'lc )  =  Prc  J  (c<<f'  I  c )  = 


(2.7) 
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=  2^X1''  ~  ^ 


d(^ 


changing  integration  order  and  integrating  with  respect  to  r  leads  to: 


WtT  J{  iL  ""P  f  T  cxp(-4)  -exp(-^-) 


(2.8) 


where: 


V  = 


(1  -c)3  , 


^r/(y)  -  er/ (ri) 


d(J) 


a: 


^  1 
P  (1-0- 


m  =  1  -csin-((j))  ,  n  =ccos-(<{)) 

y=^fW  '  Ti  =  Vr7o 


The  result  above  is  slightly  bit  different  than  that  of  Sakai  &  Tokumaru^*^. 

Differentiating  Eq.  (2.7)  with  re.spect  to  cf'  leads  to  the  probability  density  func¬ 
tion  of  the  SMC  estimator  P(("ic)  : 


rfij) 


(2.9) 


Plots  of  the  distribution  function  and  the  density  function  for  vai’ious  values  of  (f ' ,  and  the 


ratio  ~  are  enclosed.  Figuress  2.2  -  2.7  give  the  probabilty  density,  the  comulative  distri¬ 


bution,  the  mean  and  the  squared  standard  deviation  of  the  estimated  SMC  values  as  a 
function  of  the  real  SMC,  for  large  values  of  the  ratio  y  (  2048,  1024,  and  512  respec- 
tevly). 


Cam.  Prob.  Dis'tr.of  Eat.  5nrc 


fig  2.2  :  probability  density  and  comulative  distribution  of  the 
timated  SMC  as  a  function  of  the  real  SMC  for  =  2048 
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EXP.  VALUE  of  S.M.C  ESTIMATION 


VflRlRNCE  of  5.M.C  ESTIMRTION 


Figure  2.3  :  Mean  and  squared  standm’d  deviation  of  the  estimat¬ 
ed  SMC  as  a  function  of  the  real  SMC  for  y  =  2048. 
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EXP.  VPLUE  of  S.M.C  ESTIMPTION 


VARIANCE  of  S.M.C  ESTIMATION 


Figure  2.5:  Mean  and  squai'ed  standard  deviation  of  the  estimat¬ 
ed  SMC  as  a  function  of  the  real  SMC  for  =  1024. 


Com.  Prob.  Dfstr.of  Es't i.  Smc 


fig  2.6  :  probability  density  and  comulative  distribution  of  the 
timated  SMC  as  a  function  of  the  real  SMC  for  -^  =  512 
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EXP.  VALUE  of  S.M.C  ESTIMATION 


VflRlRNCE  of  S.M.C  ESTIMRTION 


Figure  2.7;  Mean  and  squared  standard  deviation  of  the  estimat¬ 
ed  SMC  as  a  function  of  the  real  SMC  for  y  =  512. 
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n.3  OTHER  TECHNIQUES  OF  SMC  ESTIMATION 


Nuttall^^  suggested  a  technique  which  enables  estimation  of  the  coherence 
directly.  For  two  given  data  sequences  xUi)  and  y(n),  two  filters  are  used  to  lineai'ly  esti¬ 
mate  y  from  x  and  also  x  from  y.  The  product  of  the  two  filter  ti'ansfer  function  is  then 
used  as  an  estimate  of  the  SMC.  In  this  way,  no  denominator  estimates  (Eq.  (1.1))  ai'e 
required. 

For  stationaiy  processes  x  and  y,  the  linear  filter  with  input  x  that  minimizes  the 
mean  square  eiTor  in  predicting  the  value  of  y  at  the  same  time  instant  has  a  transfer  func¬ 
tion  given  by  .  The  linear  predictor  of  x  based  upon  data  recorded  y  is  accom- 

pushed  by  a  linear  filter  with  mmsfer  function  .  The  product  of  these  two  transfer 

Uyy  VCO) 

function  is  then  precisely  the  SMC  between  x  and  y. 

The  linear  estimate  of  y„  from  the  available  data  record  <  x„  }  and  the  linear  esti¬ 
mate  of  x„  from  the  data  record  {  y„  }  are  given  by 

fn=±a,x„.,  (2.10) 

k=-q 

~  ^m-k 

k^p 

In  Chapter  3  and  Appendix  A  we  give  some  simulation  results  of  Nuttall's  SMC  estima¬ 
tion  method. 


n.4  THE  EFFECT  OF  SHORT  DATA  ON  SMC  ESTIMATION 


As  was  described  in  the  previous  chapters  two  main  methods  of  SMC  estimation 
have  been  used.  One  of  them  gives  the  analysis  in  the  frequency  domain,  where  the  the 
periodogram  plays  a  fundamental  role  in  the  cross  spectrum  matiix  estimation.  The  other 
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gives  the  analysis  in  the  time  domain,  by  which  we  mean  that  we  postulate  some 
parameti'ic  model  and  the  data  ai'e  fitted  to  this  model  by  estimating  the  parameters.  The 
knowledge  of  these  pai'ameters  enables  us  to  find  dh'ectly  the  SMC  function^^. 

Although  the  analysis  described  by  Carter^,  Amos  and  Koopmans^  and  Hannan*® 
holds  only  in  the  limit  [for  N  ],  for  all  practical  purposes  for  N  greater  then  40  the  law 
of  large  numbers  holds  and  it  can  be  assumed  that  this  analysis  is  valid.  It  is  obvious  that 
the  variance  of  the  estimator  and  the  bias  becomes  larger  as  the  number  of  points  decrease, 
and  hence  the  confidence  level  increases.  Estimation  of  the  SMC  in  the  frequency  domain 
is  very  sensitive  to  the  available  data  length.  From  Figure  2.1  (Caiter^^  ),  we  see  that  for  a 
relatively  small  number  of  time  segments,  the  confidence  interval  is  relatively  large. 

Brockwell  and  Davis  ^  [Chapter  10.8]  proved  that  for  a  one  dimensional  AR  pro¬ 
cess,  the  AR  model  parameters  <u'e  jointly  Gaussian  in  the  limit.  Following  their  proof,  it  is 
easy  to  show  that  their  result  can  also  be  expanded  to  the  d  dimensional  case,  i.e: 


where 


and 


([4]  -  [a })  =>  Normal 


o.duin-') 


Ro  Rq 
R-i  R-i 


ir]  = 


R-p+\  R-p+i 


Ro 

R-i 


(2.11) 


[L]  -  y?0  +  ^  R 

For  all  practical  purposes  for  N  greater  then  40,  the  law  of  large  numbers  holds, 
and  so  Eq.  (2.11)  holds  also.  Hence,  for  short  data,  the  statistics  of  the  model  estimated 
coefficients  remains  jointly  Gaussian. 

In  the  time  domain  analysis,  the  statistics  of  the  SMC  estimator  depends  only  on 
the  statistcs  of  the  model  estimated  coefficients.  As  these  statistics  rem;iin  jointly  Gaussian 
for  short  data,  the  statistics  of  the  estimated  SMC  does  not  change  also,  and  the  analysis 
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made  by  Sakai^®  holds. 

Figures  2.8  -  2.13  gives  the  density  function,  the  distribution  function,  the  squared 
standard  deviation  (  o-  )  and  mean  value  of  the  SMC  estimation  as  a  function  of  the  real 
SMC  and  the  ratio  y  (  N  is  the  data  length,  and  p  is  the  estimated  model  order  ),  for  a 
short  data  length  ( i.e  the  statistics  of  the  estimated  SMC  as  a  function  of  the  real  SMC  for 
^  =  32,16,8  ) 
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Prob  .Dens.  Func.  of  Est.Smc 


Cam.  Prob.  Ofsl:r.of  Est .  5mc 


Figure  2.8.  Probability  density  and  cumulative  distribution  of  the 

estimated  SMC  as  a  function  of  the  real  SMC  for  —  =  32 . 

P 
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EXP.  VALUE  of  S.M.C  ESTIMATION 


VflRIRNCE  of  5.M.C  ESTIMRTION 


Figure  2.9.  Mean  and  sqiuu'ed  stand;u'd  deviation  of  the  estimat¬ 
ed  SMC  as  a  function  of  the  real  SMC  for  =  32. 


Figure  2.10.  Probability  density  and  cumulative  distiibution  of 
the  estimated  SMC  as  a  function  of  the  real  SMC  for  —  =  16  . 
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EXP.  VRLUE  of  S.M.C  ESTIMATION 


VflRinNCE  of  5.M.C  ESTIMRTION 


Figure  2.11.  Mean  and  squared  standai'd  deviation  of  the  estimat¬ 
ed  SMC  as  a  function  of  the  real  SMC  for  y  =  16. 


24 


EXP.  VALUE  of  S.M.C  ESTIMATION 


VflRIRNCE  of  S.M.C  ESTIMfiTION 


Figure  2.13.  Mean  and  squai'ed  standm'd  deviation  of  the  estimat¬ 
ed  SMC  as  a  function  of  the  real  SMC  for  —  =  8. 
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m  SIMULATION  RESULTS 


Some  simulations  were  run  in  order  to  compare  the  performance  of  various 
methods  of  SMC  estimation  .  Two  models,  which  are  most  often  used  in  the  literature, 
were  chosen.  The  first  model  (denoted  as  model  1)  was  a  first  order  (p  =  1)  two  channel 
autoregressive  process  which  was  used  to  generate  data  ensembles.  Data  was  generated 
according  to  the  recursion 

+  (3.1) 

where: 


[ail  = 


0.85  -0.75 
0.65  0.55 


and  n„  is  a  zero  mean  white  Gaussian  noise  with  the  following  covariance  mauix 


Figure  3.1  gives  the  power  spectra  of  t,'  and  for  the  first  model. 


26 


Figure  3.1.  Tlic  power  spectra  of  the  signals  of  the  first  model. 


The  second  model  (denoted  as  model  2)  is  the  model  used  by  Sakai  and 
Tokumaru^®’^^ 


(3.2) 


where 


[ai]  = 


0.3  0 

0  -0.5 


,[021  = 


0.2  0 

2.0  0 


,  [03]  = 


0.1  0 
0.5  0 


and  n„  is  a  zero  mean  white  Gaussian  noise  with  the  following  covariance  matrix 


[I]  = 


1  0 
0  1 


Figure  3.2  gives  the  power  spectra  of  .t„'  and  for  the  second  model 
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Figure  3.2.  The  power  spectra  of  the  signals  of  the  second  model. 

In  the  simulation  that  we  ran,  three  methods  of  SMC  estimation  were  examined. 
The  first  was  the  traditional  way  using  the  FFT,  the  second  was  the  method  suggested  by 
Nuttall^^,  and  the  third  was  via  modeling  the  process,  estimating  its  coefficients  and  using 
Eq.  (2.5)  and  (2.6) 

The  simulation  results  are  given  in  appendix  A.  Figure  A.l  and  Figure  A.2  show 
the  true  SMC  conesponding  to  the  models  defined  by  Eq,  (3.1)  and  (3.2)  respectively.  Fig¬ 
ure  A. 3  gives  the  estimated  SMC  of  the  process  defined  by  Eq.  (3.1)(model  No.  1),  when 
the  FFT  method  was  used. 


For  the  case  when  FFT  method  was  used,  the  given  data  was  divided  into  50% 
overlapped  128  points  segments.  On  each  individual  data  segment,  the  SMC  was  calcu¬ 
lated  (using  Hamming  window).  The  estimated  SMC  is  the  average  calculated  SMC  ever 


all  data  segments.  Fora  given  data  length  of  N  points,  the  number  of  averages  (nj)  is  : 

In  order  to  have  sufficient  data  on  the  statistical  behavior  of  the  estimated  SMC,  1000 
experiments  were  nm  on  nonoverlapped  measured  data. 
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Figure  A. 4  shows  the  squared  standard  deviation  and  the  normalized  squared 
standard  deviation  of  the  estimator.  The  normalized  squai'ed  standiu'd  deviation  was 
defined  as 

var  (S -(E(<f))- 
c 

Figures  A.5  -  A. 20  give  the  histogram  of  the  estimated  SMC  at  16  different  normalized 
frequencies  co,  such  that: 

-  _  1  2  3  15  16 

-  -JI'  TI’  IT . TT:  17 

Figures  A.21  -  A.38  and  Figures  A. 39  -  A.56  are  the  same  as  Figures  A.3  -  A. 20,  but  for 
segment  length  of  256  and  128  points,  respectively. 

Figures  A.57  -  A. 188  show  the  results  of  SMC  estimation  using  Nuttall's  method, 
for  data  length  of  512,  256,  and  128  points.  Following  his  simulation  we  chose  p  =  q  =  (•> 
(Eq.  2.10).  In  order  to  get  the  statistics  of  the  estimator,  1000  Monte  Carlo  experiments 
were  run.  Figure  A.57  shows  the  results  of  all  the  runs,  for  data  segments  of  512  points 
each.  Figure  A.57.1  gives  the  averaged  estimated  SMC,  averaged  all  over  the  1000  exper¬ 
iments,  Figure  A.58  shows  the  squared  standard  deviation  and  the  normalized  squared 
standard  deviation  of  the  SMC  estimator  Figures  A. 59  -  A,74  give  the  statistics  of  the 
estimated  SMC  for  16  different  frequencies  (the  same  as  in  Figures  A.5  -  A. 20)  Figures 
A.75  -  A.87  and  Figures  A. 88  -  A. 100  give  the  statistics  of  the  two  filters  coefficients  (the 
a’s  and  the  b's  in  Eq.  (2.10))  Figures  A. 101  -  A. 144  and  Figures  A. 145  -  A. 188  iU’e  the 
same  as  Figures  A.57  -  A.  100  but  for  data  length  of  256  and  128  points  respectively. 

Figures  A.  189  -  A. 260  give  the  simulation  result  of  the  SMC  estimation  following 
Eq.  (2.6)  where  model  No.  1  was  taken.  100  nin  of  each  experiment  were  done.  Four 
cases  where  examined : 

a)  -  Data  length  of  1024  points  was  taken,  "  estimated  "  model  order  was  assumed 
to  be  1  (the  real  model  order) 
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b)  -  Data  length  of  256  points  was  taken,  "estimated"  model  order  was  assumed  to  be 

1 

c)  -  Data  length  of  128  points  was  taken,  "  estimated  "  model  order  was  assumed 

to  be  1 

d)  -  Data  length  of  512  points  was  taken,  "  estimated  "  model  order  was  assumed 

to  be  4 

The  reason  of  choosing  case ’d’  is  the  desire  to  compare  the  estimated  SMC  statistics  with 
the  result  of  case  ’b’  and  to  check  if  the  statistics  of  the  estimated  SMC  depends  only  on 
the  ratio  y  as  in  Eq.  (2.8)  and  not  in  the  values  of  N  and  p  separately.  Both  cases  'b'  and 

’d’  have  the  same  ratio  — . 

Figures  A.  189  -  A.206  show  the  average  estimated  SMC,  the  squared  standard  deviation 
the  normalized  squared  standai'd  deviation  and  the  histograms  of  the  estimated  SMC  in 
several  frequencies  for  the  case  where  1024  points  data  length  was  taken  and  model  order 
of  1  was  assumed.  Figures  A.207  -  A.224  shows  the  case  where  256  point  data  length  was 
chosen,  and  model  order  of  1  was  assumed.  Figures  A.225  -  A.227  give  the  results  for  128 
point  data  length  and  model  order  was  assumed  to  be  1.  Figure  A. 243  -  A.258  summarize 
the  case  when  data  length  of  512  points  was  taken,  and  model  order  was  assumed  to  be  4. 
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IV  DISCUSSION  ON  RESULTS 


The  statistics  of  the  model  estimated  coefficients  were  calculated  using  the  simula¬ 
tion  results.  The  coefficients  were  calculated  for  the  cases  when  method  of  estimating  the 
coherence  in  the  time  domain  was  used,  or  when  Nuttall’s  method  is  used.  The  results 
show  that  the  statistics  of  the  estimated  coefficients  remains  almost  Gaussian  even  for  rela¬ 
tively  short  data  (Figures  A.163  -  A. 188),  a  result  that  was  expected. 

Looking  at  the  probability  density  function  of  the  estimated  SMC  obtained  either 
from  the  analytic  expression  based  on  Sakai’s  work^®  [Figures  2.2  -2.7]  or  from  simulation 
results  based  on  SMC  estimation  using  the  frequency  domain  (the  FFT  method)(on  model 
No.  1),  the  time  domain  (modeling  the  process  as  an  AR  process),  or  using  NuttalFs 
method  gives  roughly  the  same  behavior  of  the  estimated  SMC  probability  density  func¬ 
tion  as  a  function  of  the  real  SMC, 

For  SMC  values  between  0.3  and  0.8,  the  probability  density  function  is  sym¬ 
metric  and  Gaussian-like.  Hence  the  assumption  of  Brockwell  and  Davis^,  and  of  Bendat 
and  Piersol^  about  the  Gaussian  behavior  of  the  density  function  when  they  calculated  the 
confidence  level  of  the  estimated  SMC  was  reasonable.  For  low  values  of  SMC  (less  then 
0.3)  and  for  high  values  of  SMC  (above  0.8)  the  probability  density  function  is  totally  dif¬ 
ferent  from  Gaussian.  It  is  an  asymmetric  function  where  the  asymmeti'y  becomes  larger 
and  larger  as  the  values  of  the  SMC  become  closer  to  zero  or  to  one.  For  these  values  of 
SMC,  BrockwelFs  or  Bendat’s  method  of  calculating  the  confidence  interval  gives  an 
inaccurate  answer. 

The  behavior  of  the  estimated  SMC  density  function  remains  roughly  the  same 
even  when  relatively  short  data  is  used,  but  as  expected,  the  standard  deviation  of  the 
estimated  coherency  became  larger.  Figure  4.1  gives  the  squared  standard  deviation  (o-) 
of  the  estimated  SMC  of  the  process  defined  by  model  No.  1  (Eq.  3.1)  when  the  method 
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based  on  the  FFT  was  used  to  estimate  the  SMC. 


SMC  Eslimaljan,  uscng  FFT,  doto  length  S12,  256  points 


Figure  4.1.  Valuation  of  squared  standard  deviation  of  estimated 

SMC  as  a  function  of  data  length  (model  No.  1,  using  FFT  method). 

When  the  FFT  method  is  used  to  estimate  the  SMC,  it  can  be  seen  that  the  estimator  is 
biased  and  the  bias  increases  as  the  data  length  decreases.  This  phenomena  is  strongly 
dominant  for  low  values  of  coherence  where  the  bias  eiTor  is  much  larger  then  the  real 
value  of  the  SMC  itself.  Moreover,  caieful  inspection  in  Figure  4.2  shows  that  the  bias  (as 
well  the  standai'd  deviation)  does  not  depend  only  on  the  real  SMC  value  but  also  on  the 
frequency,  a  result  which  is  not  given  in  the  expression  by  Cai'ter  and  Sakai. 
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SMC  esilmation,  using  FFT 


Figure  4.2.  Bias  vai  iation  of  the  estimated  SMC  as  a  function  of  data  length 
(model  No.  1,  using  FFT  method) 


Figure  4.3  gives  the  simulation  results  and  calculated  results  of  the  estimated  SMC 
variance  for  data  length  of  256  points  (i.e  nj  =  3).  The  calculated  result  was  obtained  by 
using  the  approximation  given  by  Carter  et  al^.  As  expected,  in  the  region  where  the 
statistics  of  the  estimated  SMC  aie  Gaussian-like,  the  calculated  and  the  simulated  vari¬ 
ance  are  relatively  similar.  For  low  values  of  SMC,  the  statistics  of  the  estimated  SMC  is 
not  Gaussian  any  longer  and  the  calculated  and  the  simulated  vai’iance  are  different.  As  a 
conclusion,  we  see  that  the  calculated  variance  is  a  relatively  good  upper  bound  on  the 
estimated  SMC  vai’iance  for  SMC  values  ku-ger  then  0.3  even  for  very  short  data. 
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SMC  esiltnoilon,  <-<5  tng  FFT,  segment  length  256  points 


Figure  4.3.  Compaiison  of  the  calculated  and  simulated  o-  of  the 

estimated  SMC,  for  nrf=3,when  FFT  method  is 
•  -■  used. 

For  model  No  1,  assuming  model  order  as  the  real  order  (1),  a  compiuison  of  the 
SMC  estimation  results  between  the  "FFT"  method,  NuttalFs  method  and  the  "direct" 
method  based  on  the  frequency  domain,  was  made.  The  compiuison  shows  that  the  direct 
method  (Eq.  2.4  -  Eq.  2.6)  is  superior  in  the  sen.se  of  bias  eri'or  as  well  as  the  vaiaance  of 
the  estimation.  The  worse  method  is  the  FFT  method  which  has  the  largest  bias  and  viui- 
ance.  The  difference  in  the  performance  between  the  three  methods  is  emphasized  espcci- 
aly  at  low  values  of  SMC.  Figure  4.4  shows  the  average  SMC  estimation  by  the  three 
method  for  model  No  1,  model  order  taken  as  1  (the  real  order),  and  data  length  of  128 
samples. 
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SMC  estimolian,  d^to  length  I  ?B  points 


Figure  4.4.  Comparison  of  die  average  estimated  SMC  for  model  No  1, 

data  lengtli  128  samples,  and  assumed  model  order 
is  1 . 


Examing  the  results  of  the  SMC  estimation  when  model  No  2  is  used  and  compar¬ 
ing  it  with  the  results  when  model  No  1  was  taken  leads  to  the  same  conclusion  as  before: 
the  direct  method  is  the  best  and  the  FFT  is  the  worst.  An  interesting  point  is  that  the 
values  of  the  bias  and  the  variance  were  slightly  different,  even  though  same  data  length 
was  taken,  or  different  data  lengths  were  taken  but  such  that  the  ratio  y  was  the  same.  The 
comparison  showes  that  the  SMC  estimation  is  worst  for  model  No  2. 

One  of  the  results  we  obtained  following  Sakai  and  Tokumaru^^  is  that  the  statis¬ 
tics  of  the  estimated  SMC  depend  among  other  parameters  on  the  ratio  between  the  data 
length,  N,  and  the  model  order,  p.  Comparing  the  simulation  results  of  two  cases,  'b'  (Fig¬ 
ures  A. 207  -  A.225)  and  'd'(  Figures  A.243  -  A. 270  )  which  have  the  same  ratio  ^ 


35 


(y  =  128),  leads  to  a  different  conclution.  The  statistics  of  the  estimated  SMC,  the  variance 

and  the  bias  in  these  two  cases,  are  slightly  different  (compare  Figure  A. 207.1  with 
A.243.1  and  Figure  A.208.1  with  A. 244.1  ),  which  leads  to  the  conclusion  that  the  statis¬ 
tics  of  the  estimated  SMC  depend  not  only  on  the  ratio  y  but  also  on  the  parameters  N 
and  p  themselves. 
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APPENDIX  A 


flvorege  S.M.C  (using  FFT) 


I 
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Normalized  Frequency 

Fjg.  fi.3  ;  SMC  es-timaiion,  using  FFT  method,  data  length  512  points, 
Rverage  estimated  SMC,  over  1000  runs 


Nonrolized  5igiia*<2  of  S.M.C  (using  FR) 


0.0  0.1  0.2  0.3  0.4  0.5 

Normalized  Frequency 

SMC  estimaiion,  using  FFt  meihod,  data  length  512  points, 

Squared  Standard  Deviation  of  estimated  SMC  data  length  512  points,  1000 


,  UMM  rxr«. 


Normalized  Frequency 

Ei3.-__Q.-A-_1_!  SMC  estimation,  using  FFT  method,  data  length  512  points, 

Normolized  squared  stancfcrd  deviation  of  estimated  SMC  over  1000  runs 


rlisTogram  of  Coherency  (via  FFT)  at  freq. 
Mean=0. 31533, Var=0. 33912e -01 


=  1/32 


0.09  0.17  0.25  0.35  0.43  0.52  0.61  0.69  0.78  0.87 

Bin  Centers 

Ei9.-_Jl-_5_u  SMC  estimaiion,  using  FFT  method,  data  length  512  points, 
Histogram  of  estimated  SMC  ot  normalized  frequency  =  1/32 


Histogram  of  Coherency  (via  FFT)  at  freq. 
Mean=0 .61198, Var=0 . 28391e-01 


=  2/32 


0.09  0.19  0.26  0.37  0.47  0.56  0.66  0.75  0.64  0.94 

Bln  Cen+ers 

SMC  estimation,  using  FFT  method,  data  length  512  points. 
Histogram  of  estimated  SMC  at  normalized  frequency  =  2/32 


3/32 


Cbserwtions/Bin  Cbsavotiots/Bin 


Hls-togram  of  Coherency  (vlo  FFT)  of  freq.  =  5/32 


Mean=0. 96151 , Var=0. 23971b-01 
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Bin  Cenier-s 

SMC  estimaiion,  using  FFT  method,  data  length  512  points 
Histogram  of  estimated  SMC  ot  normalized  frequency  =  5/32 


Histogram  of  Coherency  (via  FFT)  at  freq. 
Mean=0. 64600, Var=0. 66576e-01 


=  6/32 


0.41  0.47  0.53  0.59  0.65  0.71  0.77  0.83  0.89  0.95 

Bln  Ceniers 

Fi9^.__fl._l_0_!  SMC  estimation,  using  FFT  method,  data  length  512  points 
Histogram  of  estimated  SMC  at  normalized  frequency  =  6/32 


7/32 


rllSTogram  of  Cohonency  (via  FFT)  of  freq 
Mean=0 . 47135 . Var=0 . 19310 


=  9/32 
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Fig.  P.13  ;  SMC  estimation,  using  FFT  method,  data  length  512  points. 
Histogram  of  estimated  SMC  at  normalized  frequency  =  9/32 


Hi s-togrom  of  Coherency  (via  FFT)  at  freq. 
Mean=0. 33251 , Var=0. 16719 


=  10/32 
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Fig.  R.21  !  SMC  estimo-tion,  using  FFT  method,  data  length  256  points, 
Average  estimated  SMC  over  1000  runs 
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Eia--  jL._22_A:  SMC  estimo' 
Normalized  squoi 


1/32 


His-togram  of  Coherency  (vfo  FFT)  ot  freq.  = 
Mean=0 . 440&9 , Var=0 . 65&10a-01 
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Fig.  P.23  ;  SMC  estimation,  using  FFT  method,  data  length  256  points. 
Histogram  of  estimated  SMC  at  normalized  frequency  =  1/32 


Histiogrom  of  Coherency  <vio  FFT)  of  freq.  =  2/32 
Mean=0. 67301 . Var=0. 45800e-01 
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Bln  Centers 

Fig..__fi._24 SMC  estimation,  using  FFT  method,  data  length  256  points, 
Histogram  of  estimated  SMC  at  normalized  frequency  =  2/32 


3/32 


Histogram  of  Cohsrsncy  (via  FFT)  at  freq.  = 
Mean=0. 88903. Var=0. 12180e-01 


0.18  0.27  0.36  0.45  0.54  0.63  0.72  0.81  0.89  0.98 

Bin  Cen-ter-s 

Eia._p  ,_25_!  SMC  estimation,  using  FFT  method,  data  length  256  points, 
Histogram  of  estimated  SMC  at  normalized  frequency  =  3/32 


Histogram  of  Coherency  (via  FFT)  at  freq.  =  4/32 
Mean=0 . 99725 . Var=0 . 14475e-04 
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Bln  Centers 

Fia.--a  ._26_!  SMC  estimation,  using  FFT  method,  data  length  256  points, 
Histogram  of  estimated  SMC  at  normalized  frequency  =  4/32 


5/32 


n.i8 


Histogram  of  Cohsrsncy  (via  FFT)  at  frsq.  =  7/32 

MeQn=0 . 7A1S^ ,  Var=0 . 392E53e-01 


Bln  Ceniers 

EiQ.-__fi-_29 SMC  estimation,  using  FFT  method,  data  length  256  points. 
Histogram  of  estimated  SMC  ot  normalized  frequency  =  7/32 


Histogram  of  Coherency  (via  FFT)  at  freq.  =  8/32 
Mean=0. 65741 , Var=0. 54133e-01 


Bln  Cen+ei-s 

Fi3..__R._30_!  SMC  estimation,  using  FFT  method,  data  length  256  points. 
Histogram  of  estimated  SMC  at  normalized  frequency  =  8/32 
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Histogr-am  of  Coherency  (via  FFT)  oi  freq.  =  9/32 


Mean=0 . 5"7301  ,  Var'=0 . 63562b-01 
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Fig.  R.31  !  SMC  estimation,  using  FFT  method,  data  length  256  points. 
Histogram  of  estimated  SMC  at  normalized  frequency  =  9/32 


Hisiogrom  of  Coherency  (via  FFT)  o-t  freq 
Mean=0. 51373, Var=0. 64967e-01 


=  10/32 
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Fig..__0.._32 SMC  estimation,  using  FFT  method,  data  length  256  points, 
Histogram  of  estimated  SMC  at  normalized  frequency  =  10/3i 
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Fig^.__fl._40^._l:  SMC  estimotion,  using  FFT  method, 
Normalized  squared  standard  deviation  c 
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His-togram  of  Coherency  (via  FFT)  at  freq.  =  1/32 
Mean=0 . 40700 , Var=0 . 63370e-01 


Bln  Ceniers 

Ei3.._Jl-_4L_!  SMC  estimation,  using  FFT  method,  doto  length  128  points. 
Histogram  of  estimated  SMC  ot  normalized  frequency  =  1/32 


Histogram  of  Coherency  (via  FFT)  at  freq.  =  2/32 
Mean=0. 61 133. Var=0. 56663e-01 


Bln  Ceniers 

Fig_.__fi._42_!  SMC  estimotion,  using  FFT  method,  data  length  128  points, 
Histogram  of  estimoted  SMC  at  normalized  frequency  =  2/32 


Cbservotions/Bin  Obsendions/Bin 


HlsTognam  of  Coherency  (vio  FFT)  of  freq. 
Mean=0 . 64649, Var=0.  16714b-01 


=  3/32 
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Fig.  R.43  !  SMC  estimation,  using  FFT  method,  data  length  128  points, 
Histogram  of  estimated  SMC  at  normalized  frequency  =  3/32 


Histogram  of  Coherency  (via  FFT)  at  freq. 
Mean=0. 99246, Var=0. 94314e-04 


=  4/32 
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Bin  Ceniens 

F{g^.__fi._44__:  SMC  estimation,  using  FFT  method,  data  length  128  points 
Histogram  of  estimated  SMC  at  normalized  frequency  =  4/32 


5/32 


Histogram  of  Coherency  (via  FFT)  at  freq. 
Mean=0. 77360, Var=0.  34611a-01 


=  7/32 
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Fig.  P.47  ;  SMC  estirnDiion,  using  FFT  method,  data  length  128  points 
Histogram  of  estimated  SMC  ot  normalized  frequency  =  7/32 


Histogram  of  Coherency  (via  FFT)  at  freq. 
Mean=0 . 67161 , Var=0 . 50669e-01 


=  8/32 
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Fig..__fi._48 SMC  estimation,  using  FFT  method,  data  length  128  points, 
Histogram  of  estimated  SMC  at  normalized  frequency  =  8/32 
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ObservnIiaRs/Bin  Observntions/Bin 


His+ogram  of  Coherency  ox  normalized  f  req .  =:  11/32 

MeQn=0 . 22773 ,  Squared  Standard  Deviatlon=0.  66822e-02 
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BJn  Centers 

El3.-__0.-_69 _!  SMC  estimation,  Nuttall’s  method,  data  length  512  points, 
Histogram  of  estimated  SMC  at  normalized  frequency  =  11 /3c 


Histogram  of  Coherency  at  normalized  freq.  =  12/32 
Mean=0 .  ■*  6‘^69 ,  Squared  Standard  Devlot Jon=0. 56616e-02 


0.05  0.10  0.16  0.21  0.26  0.31  0.36  0.41  0.47  0.52 

Bln  Centers 

Ei9.-__0.-_Z9  SMC  Gstimotion,  Nuttall's  method,  data  length  512  poin 
Histcgrom  of  estimoted  SMC  at  normalized  frequercy  =  12 


Hisiognam  of  Coherency  of  normalized  freq.  =  13/32 
Mecin=0 .  1 1 230 ,  Squared  S'tandard  Devlot Jon=0.  35522e-02 


0.04  0.07  0.11  0.15  0.18  0.22  0.26  0.29  0.33  0.37 

Bin  Cen"ker& 

Ei3.-__Q.JZ.l-_!  SMC  estimation,  Nuttall’s  method,  data  length  512  points, 
Histogram  of  estimated  SMC  at  normalized  frequency  =  13/3£ 


Hisfogrom  of  Coherency  of  normalized  freq.  =  14/32 
Mean=0 . 66639e-01 ,  Squared  Sfondord  Devia‘txon=0. 25034e-0, 
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Ei3.-__Q.-_Z2 _:  SMC  estimation,  Nuttall’s  method,  data  length  512  points 
Histogram  of  estimoted  SMC  at  normalized  frequency  =  14/3 


His-tognam  of  Coherency  of  normolized  freq.  =  15/32 
Mean=0 . 45600b -01  ,  Squared  Stondard  Devia'tfon=0.  15243b-02 
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Ei9.-_J0.J2L_!  SMC  estimation,  Nuttall's  method,  data  length  512  points, 
Histogram  of  estimated  SMC  at  normalized  frequency  =  15/32 

Hisfogrom  of  Coherency  of  normalized  freq.  =  16/32 
Mean=0. 382506-01 ,  Squared  Standard  Devlation=0. 21616e-02 


0.03  0.07  0.10  0.13  0.17  0.20  0.24  0.27  0.30  0.34 


Bin  Centers 

Fi^.__R._74_:  SMC  estimation,  Nuttall’s  method,  data  length  512  points. 
Histogram  of  estimated  SMC  at  normalized  frequency  =  16/3^ 
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Mean=-0.  1  497’6e-01  .  Siandard  Devfa‘tIon=0.  34077e-01 
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SMC  estimation,  Nuttall's  method,  data  length  512  points. 
Histogram  of  estimated  model  coeficient  ’P)4’ 


His-togram  of  R4 

Mean=-0.  44590B-01  ,  Standard  DevfDtIon=0.  34920e-01 


Bin  Centers 

_=  SMC  estimation,  Nuttoll’s  method,  data  length  512  points, 
Histogram  of  estimated  model  coeficient  '03' 
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Fig,  fl.65  i  SMC  estimation,  Nuttoll's  method,  data  length  512  points, 
Histogrom  of  estimated  model  coefjcient 


Histogram  of  fll2 

Mean=0 . 55921e-02 ,  Standard  DBVlatIon=0.  32301b-01 
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Fig^.__fi._86 _!  SMC  estimotion,  Nuttoll's  method,  data  length  512  points, 
Histogram  of  estimated  model  coeficient  ’fill 


Cbservolians/Bin  CbsBMitjons/Bin 


Cbservotions/Bin  ObsrwtJore/Bin 


Histogram  of  B5 


Mean=0. 22547.  siandard  devJatJon=0. 33V70e-01 
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SMC  estimation,  Nuttall’s  method,  data  length  512  points. 
Histogram  of  estimated  model  coeficient  ’B6’ 


Histogram  of  B6 

Mean=0 . 46214 ,  standard  devlat  Jon=0 . 3921  le  -  01 
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Fig^.__fi._93 _t  SMC  estimation,  Nuttoll’s  method,  data  length  512  points, 
Histogram  of  estimated  model  coeficient  ’B5’ 
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Mean=0  .  17790 ,  standard  devlatlon=0 . 36951e“01 
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Fi£. _ fi._94_j  SMC  estimotjon,  Nuttoll’s  method,  data  length  512  points, 

Histogram  of  estimated  model  coeficient  ’Bfl’ 


Hisfogram  of  B8 

Mean=-0. 35241  ,  standard  devlat/Dn=0.  37902e-01 
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Fl2^.__R._95  SMC  estimotion,  Nuttoll’s  method,  doto  length  512  points 
Histogram  of  est/moted  model  coeficient  ’B7’ 
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Hls-fcognam  of  B13 

Mean=-0.  15034e-01  ,  siandar'd  de\/iation=0 . 26165e-01 
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Fj'g.  R.  100  !  SMC  estimation,  Nuttall’s  method,  data  length  512  points. 
Histogram  of  estimated  model  coeficient 
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HlsTor-am  of  Coherency  ot  normalized  freq.  =  11/32 

Mean=0 . 2435S ,  Siandard  Devlatfon=0.  12626e-01 
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:  SMC  esiimaiion,  Nutioll’s  method,  data  length  256  points, 
Histogram  of  estimated  SMC  at  normoliiied  frequency  =  11 /3c 


Histogram  of  Coherency  ot  normolized  freq. 

Mean=0.  18550 ,  Standard  Devlatlon=0.  10737e-01 


=  12/32 
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Ei3.-__8.-_l.L4  !  SMC  estimation,  Nuttall’s  method,  data  length  256  points 
Histogrom  of  estimated  SMC  at  normalized  frequency  =  12/C 
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Histogram  of  R3 

Mean=-0.  1  5421b- 01  ,  S'tandard  Dovfa-tlon=0.  45&'72e>-01 


-0.11-0.08-0.06-0.03  0.00  0.03  0.06  0.08  0.11  0.14 

Bln  Cen-lers 

i  SMC  estimation,  Nuttoll’s  method,  data  length  256  points 
Histogram  of  estimated  model  coeficient  'fi4' 


His-togrom  of  R4 

Mean=-0, 45569b-01  ,  Standard  Dev£at/on=0,  46792e-01 
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Ei9.-__8.-_l.?.?.  =  SMC  estimation,  Nuttall's  method,  data  length  256  points 
Histogram  of  estimated  model  coeficient  ’fl3’ 
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Fig.  P.123  I  SMC  esfimation,  Nuitall’s  method,  data  length  256  points, 
Histogram  of  estimated  model  coeficient  ’ft6’ 
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Mean=-0. 36'799,  Standard  Devlation=0. 59665e-01 
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Ei9.-_Ji-_l.24,  !  SMC  estimation,  Nuttall’s  method,  data  length  256  points. 
Histogram  of  estimated  model  coeficient  ’fl5' 
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His+ognam  of  09 

Mean=0.  14539 ,  S-tandard  DBviatIon=0. 50574e-01 
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Fi9..__fl._127  :  SMC  estimation,  Nuttall’s  method,  data  length  256  points 
Histogram  of  estimated  model  coeficient 


Histiogrnm  of  010 

Mean=0. 49067e-01 ,  Stondard  DevlatIon=0. 48060e-01 


0.15-0.11-0.06-0.04  0.00  0.04  0.08  0.11  0.15  0.19 

Bln  Cenier-s 

Fia-__a--  _^28  !  SMC  estimation,  Nuttall’s  method,  data  length  256  points 
Histogram  of  estimated  model  coeficient  ’R9’ 


Cbservoticns/Bin  Cbsavitions/Bin 


nlsTogram  of  flll 

Mean=0.  1 6069e- 01  ,  Standard  Devlation=0. 47125e-01 


-0.13  0.10-0.07-0.03  0.00  0.03  0.07  0.10  0.13  0.16 

Bln  Centers 

Fig.  P.129  s  SMC  estimation,  Nuttall’s  method,  data  length  256  points 
Histogrom  of  estimated  model  coeficient 


Histogram  of  R12 

Mean=0. 50536e-02,  Standard  DBvlatlon=0. 44996e-01 


-0.16-0.13  0.10-0.07-0.03  0.00  0.03  0.07  0.10  0.13 

Bln  Centers 

Fl£.__fl.JJ30  !  SMC  estimation,  Nuttall’s  method,  data  length  256  points 
Histogram  of  estimated  model  coeficient  ’fllr 


Hls"togi~am  of  P13 

Mean=0. 33855e-02,  Standard  Devfatlon=0. 3547’4e-01 


-0.09-0.07-0.05-0.02  0.00  0.02  0.05  0.07  0.09  0.12 


Bin  Cen-ters 

Fig.  R.  131  !  SMC  esfimaiion,  Nuttall’s  mgihod,  dcrta  length  256  points, 
Histogram  of  estimated  model  coeficient 


Cbservoticns/Bin  Cbservotions/Bin 


Hlstognam  of  B3 

Mean=0 . 464&le-01 ,  s'tandord  devla-tion=0 . 41551e-01 


-0.05-0.03  0.00  0.03  0.05  0.08  0.11  0.13  0.16  0.19 

Bln  Centera 

Fig,  ft.  134  !  SMC  esiimation,  NutioH’s  method,  data  length  256  points, 
Histogram  of  estimated  model  coeficient  ’B4' 


Histogram  of  B4 

Mean=0.  10704,  standard  devJation=0. 44966e-01 


0.00  0.03  0.06  0.09  0.11  0.14  0.17  0.20  0.23  0.26 

Bin  Cen-ters 

Ei3.-__8.JL35  :  SMC  estimation,  NuttoH's  method,  data  length  256  points. 
Histogram  of  estimated  model  coeficient  ’B3’ 


Cbsemticns/Bin  Cbsavitions/Bin 


Histogram  of  BV 

Mean=0.  17643 ,  siandard  devlot  JDn=0 . 53V15e-01 


0.03  0.07  0.10  0.14  0.17  0.21  0.24  0.27  0.31  0.34 

BJn  Cenlei-s 

Fig^.__fl._^38  !  SMC  esiimaiion,  Nutiall’s  method,  data  length  256  points, 

Histogram  of  estimated  model  coeficient  ’B8’ 


Histogram  of  BS 

Mean=-0. 35516,  standard  devJotlon=0, 56201b-01 


0.54  0.50-0. 46-0. 43-0. 39-0. 36-0.32-0.29-0.25-0.21 

Bjn  Centers 

Fig^.__fl._lJ39  !  SMC  estimation,  Nuttall’s  method,  data  length  256  points. 
Histogram  of  estimated  model  coeficient  ’B7’ 


CbserMitions/Bin  CheervotJons/Bin 


fl.82 

SMC  esi  Jmailon ,  Nuttal’s  method,  data  length  128  points 


Normalized  Frequency 

Fig.  R.146  :  Squared  Standard  Deviation,  of  estimated  SMC,  Nuttall’s  method 

data  length  128  points,  1000  runs 

SMC  estimation,  Nuttol's  method,  data  length  128  points 


Normalized  Frequency 


Fig.  fl.  146.1:  Normalized  Squored  Stondard  Deviation  of  estimated  SMC,  Nuttsll’s  method 

data  length  128  points,  1000  runs 


Histogram  of  Cohenency  oi  nor-moliZQd  froq. 
Mean=0, 23572 ,  Standard  DavIa-tlon=0.  12243 


=  1/32 


0.09  0.17  0.26  0.35  0.43  0.52  0.60  0.69  0.76  0.86 

Bin  Centers 

El9.-_J1-_1.4Z  s  SMC  estimation,  Nuttall’s  method,  data  length  128  points. 
Histogram  of  estimated  SMC  at  normalized  frequency  =  1/32 


His+ogram  of  Coherency  at  normalized  freq. 
Mean=0. 57273,  Standard  DevIatlon=0.  10629 


=  2/32 


0.21  0,26  0,35  0.42  0.49  0.56  0.63  0.69  0.76  0.63 

Bin  Centers 

Ei9.-__8.-_l.48  !  SMC  estimation,  Nuttall’s  method,  data  length  128  points. 
Histogram  of  estimated  SMC  at  normalized  frequency  =  2/32 


His-togram  of  CoherGncy  oi  normalized  freq. 
Meon=0 . 66556 ,  Standard  Devla'tlon=0. 72996e-01 


=  3/32 


0.59  0.64  0.69  0.74  0.79  0.84  0.89  0.94  0.99  1.04 

Bin  Cen-ters 

Ei9.-_J1JA9.  !  SMC  esiimation,  Nutiall’s  method,  data  length  128  points, 
Histogram  of  estimated  SMC  at  normalized  frequency  =  3/32 


His-togram  of  Coherency  oi  normalized  freq. 
Mean=0 . 99766 ,  Siandar d  Devlatlon=0.  11664e-01 


=  4/32 


0.65  0.67  0.69  0.91  0.93  0.95  0.97  0.99  1.01  1.03 

Bln  Cerrters 

Ei3.-__8.-J5.0  !  SMC  estimation,  Nuttall’s  method,  data  length  128  points, 
Histogram  of  estimated  SMC  at  normalized  frequency  =  4/32 


5/32 


Hisiognam  of  Coherency  of  normalized  freq.  = 
Mean=0 . 95420,  S-tandard  DavfaiIon=0.  28563e-01 


— 

- 

0.71  0.74  0.78  0.81  0.84  0.88  0.91  0.95  0.98  1.01 

Bln  Csnier-s 


:  SMC  estimation,  Nuttall’s  method,  data  length  128  points. 
Histogram  of  estimated  SMC  at  normalized  frequency  =  5/32 


His-logram  of  Coherency  at  normalized  freq.  =  6/32 
Mean=0. 84846,  S'tondard  Devlatlon=0.  66012e-01 


Bln  Cenier-s 


Fi^.__fl._l_52  s  SMC  estimation,  Nuttall’s  method,  data  length  128  points, 
Histogram  of  estimated  SMC  at  normalized  frequency  =  6/32 


Histogram  of  Cohensncy  oi  normalized  freq, 
Mean=0. 72012,  5-tandard  Devlatlon=0. 96322e-01 


=  7/32 


0,41  0.47  0.53  0.59  0.65  0.71  0.77  0.83  0.89  0.95 

Bin  Centers 

Ei3.-__Q.J:53.  s  SMC  estima-lion,  Nut-tail’s  method,  dato  length  128  points, 
Histogram  of  estimated  SMC  at  normalized  frequency  =  7/32 


Histogram  of  Coherency  at  normalized  freq. 
Mean=0 . 57650 ,  S-tandard  Devlatlon=0.  13149 


=  8/3^ 


0.17  0,25  0.33  0.42  0.50  0.58  0.67  0.75  0,83  0.92 

Bln  Cenier-s 

Fig..__fi._l_54  !  SMC  estimation,  Nuttall’s  method,  data  length  128  points. 
Histogram  of  estimated  SMC  at  normalized  frequency  =  8/32 


His-togram  of  Coherency  oi  normalized  freq, 
Mean=0 . 43416 ,  Siandard  DevIa'tlon=0.  14736 


=  9/32 


0.09  0.19  0.28  0.38  0.47  0.57  0.66  0.76  0.85  0.95 

Bin  Cen-ters 

Ei3.-__0.-_l.55  s  SMC  estimation,  Nuttall’s  method,  dott’  length  128  points. 
Histogram  of  estimated  SMC  at  normalized  frequency  =  9/32 


Hisiogrom  of  Coherency  ot  normalized  freq. 
Mean=0 . 33645 ,  Standard  DevfatIon=0.  14408 


=  10/32 


0.08  0.16  0.24  0.32  0.40  0.46  0.56  0.64  0.72  0.60 

Bin  Centers 

Ei3.-__0.-_l.56  :  SMC  estimation,  Nuttall’s  method,  dota  length  128  points, 
Histogram  of  estimoted  SMC  at  normalized  frequency  =  10/32 


11/32 


Hisiogram  of  Coherency  oi  normolized  freq.  = 
Mean=0 . 2'?729 ,  Standard  Devlation=0.  15153 


BJn  Centers 

Ei3.-_il-_15Z  s  SMC  estimation,  Nuttall’s  method,  data  length  128  points, 
Histogram  of  estimated  SMC  at  normalized  frequency  =  11/32 


His-togram  of  Coherency  o-t  normalized  freq.  =  12/32 
Mean=0 . 22 1 25 ,  Standard  Dev  I  at  lon=0 .  1 3793 


0.07  0.13  0.20  0.26  0.33  0.40  0.46  0.53  0.59  0.66 

Bin  Centers 


Ei9.-__Q.-_158  !  SMC  estimation,  Nuttall’s  method,  data  length  128  points. 
Histogram  of  estimated  SMC  at  normalized  frequency  =  12/32 


Hlsiogr-am  of  Coherency  oi  normolized  freq. 
Mean=0  .  17516 ,  Standard  DevlatIon=0.  11943 


=  13/32 


0.06  0.13  0.19  0.26  0.32  0.39  0.45  0.51  0.58  0.64 

Bln  Centers 

Ei3.-_J1._1.59  !  SMC  eslimation,  Nui-talTs  method,  data  length  128  pofnt'^, 
Histogram  of  estimated  SMC  at  normalized  frequency  =  13/32 


Hist:ogram  of  Coherency  ot  normolized  freq. 
Mean=0.  13714 ,  Standard  DevIatlan=0.  10666 


=  14/32 


0.06  0.11  0.17  0.22  0.26  0.33  0.39  0.45  0.50  0.56 

Bln  Centers 

Ei3.-__B.-_l.5.0  ■  SMC  estimation,  Nuttall's  method,  data  length  128  points, 
Histogrom  of  estimated  SMC  at  normalized  frequency  =  14/32 


Cbservatians/Bin  ObserwitiatB/Bin 


Histogram  of  CohsrQncy  oi  normalized  freq.  =  15/32 
Mean=0  ,  11691  .Standard  DavlatIon=0.  10475 


Bln  Centers 


i  SMC  estimation,  NuttalTs  method,  data  length  128  points. 
Histogram  of  estimated  SMC  at  normalized  frequency  =  15/32 


His-togrom  of  Coherency  at  normalized  freq.  =  16/32 
Mean=0.  11 107,  Standard  DevlatIon=0.  14262 


Bln  Centers 

Fi3..__fl._l_62  !  SMC  estimation,  Nuttoll’s  method,  data  length  128  points, 
Histogram  of  estimated  SMC  at  normalized  frequency  =  16/32 


Cbservotiona/Bin  Obaavitiore/Bin 


Hisiogram  of  h3 

Mean=-0.  1  01  ,  Standard  Devfatlon=0.  63567e-01 


-0.19-0.15-0.11-0.07-0.04  0.00  0.04  0.07  0.11  0.15 

Bin  Centers 

Fig,  ft.  165  !  SMC  estimation,  Nuttall’s  method,  data  length  128  points 
Histogram  of  estimated  model  coeficient  ’fi4’ 


His-togrom  of  R4 

Mean=-0.  49277e-01 ,  Standard  Devlatlon=0.  67606e-01 


-0.22-0. lB-0. 13-0.09-0.04  0.00  0.04  0.09  0.13  0.18 

Bin  Centers 

Fig_.__fi._1^66  !  SMC  estimation,  Nuttoll’s  method,  data  length  128  points, 
Histogram  of  estimated  model  coeficient  ’fl3’ 


Cbservotions/Bin  Cbsffvotiore/Bin 


rlisTogram  of  Ho 

Mean=-0.  13032 ,  Standard  Dev!atlon=0 . 69405e-01 


-0.32-0.28-0.23-0.19-0.14-0.09-0.05  0.00  0.05  0.09 

Bln  Centers 

••  esilmation,  Nuttall's  method,  data  length  128  points, 
Histogram  of  estimated  model  coeficient  ’fi6’ 


Histogram  of  PS 

Mean=-0. 38904,  Standard  Devlot  Jon=0. 69156e-01 


-0.67-0.61-0.55-0.49-0.43-0.37  0.30-0.24-0. 16-0.12 

Bln  Centers 

!  SMC  estimation,  Nuttall’s  method,  data  length  128  points 
Histogram  of  estimated  model  coeHcient  ’fi5' 


Cbsendiore/Bin  Cbsavilions/Sin 


rlisTogram  of  R9 

Mean=0  .  14539 ,  Standard  DevIatlon=0.  70574e-01 


0.10-0.05  0,00  0.05  0.10  0.15  0.20  0.24  0.29  0.34 

Bin  Centers 

*  SMC  esiimation,  Nuttoll’s  method,  data  length  128  points 
Histogram  of  estimated  model  coeficient 


Histogram  of  R10 

Mean=0. 46110B-01 ,  Standard  Devlatfon=0. 6©515e-01 


0.17-0,13-0.09-0.04  0.00  0.04  0.09  0.13  0.17  0.21 

Bin  Centers 

Fi3^.__R._172,  :  SMC  estimation,  Nuttall’s  method,  data  length  128  points. 
Histogram  of  estimated  model  coeficient  ’fi9’ 


Cbservoficns/Bin  Cbservjtjons/Bin 


Hls+ogram  of  B9 

Msan=-0.  17360,  standard  devlatlor>=0,  66976b-01 


-0.38-0.34  0.30-0.26-0.21-0.17-0.13-0.09-0.04  0.00 

Bln  Centers 

'  SMC  esiimation,  Nuttall’s  method,  data  lenqth  128  points, 
Histogram  of  estimated  model  coeficient  ’B10^ 


Histogram  of  B10 

Mean=-0. 84265e-01 .  standard  devlat.ion=0. 66344b-01 


0.31-0.25-0,19-0,12-0.06  0.00  0.06  0.12  0.19  0.25 

Bln  Centers 

■  SMC  estimation,  Nuttall’s  method,  data  length  128  points, 
Histogram  of  estimated  model  coeficient  ’B9’ 


Normalized  Frequency 


El9.-_J1-_189.  :  SMC  estimation,  assumed  model  order  =  1, 
data  length  1024  points 


SMC  esiimation,  general  meihod ,  data  lengih  1024  points 


0.0  0.1  0.2  0.3  0.4  0.5 

Normalized  Frequency 

Fig..__R._169. 1  :  SMC  estimotion,  assumed  model  order  =  1, 

data  length  1024  points,  flveroge  estimated  SMC  over  100  runs 


3.56  0 
B 


iima-tion 
am  of  e£ 


ilil 


Cbservoticns/Bin  Cbservotiore/Bin 


Hls-tognam  of  Coherency  at  normalized  frequecy  =  3/32 
Mean=0 . 67603  ,  Var=0 . 20951e-04 ,  assumed  model  order  =  1 


Bln  Centers 

Fig.  P.193  :  SMC  eslfmation,  assumed  model  order  =  1,  (model  1),  data  length  1024 
Histogram  of  estimated  SMC  at  normalized  frequency  =  3/32 


Histogram  of  Coherency  at  normalized  frequecy  =  4/32 
Mean=0 . 99616  ,  Var=0 . 69094e-06,  assumed  model  order  =  1 


Bln  Cen-ters 

Fig_.__fi._r94  !  SMC  estimation,  assumed  model  order  =  1,  (model  1),  data  length  102/ 
Histogram  of  estimated  SMC  at  normalized  frequency  =  4/32 


Cbservotions/Bin  Observotiora/Bin 


fl.l09 


Histogram  of  Coherency  at  normalized  frequecy  =  5/3'.' 
Mean=0 . 95622  ,  Vor=0 . 77996s -05,  assumed  model  order  =  1 


Bjn  Cerrters 

Fl£.__R._l_95  I  SMC  es-timation,  assumed  model  order  =  1,  (model  1),  doia  length  1024  polr,;: 
Histogrom  of  estimated  SMC  at  normalized  frequency  =  5/32 


Histogram  of  Coherency  at  normalized  frequecy  =  6/32 
MeQn=0. 63740  ,  Var=0 . 23607e-04,  assumed  model  order  =  1 


Bln  Cen-ters 

Fi£.__fl._1^96  :  SMC  estimation,  ossumed  model  order  =  1,  (model  1),  data  length  1024  point' 
Histogram  of  estimated  SMC  at  normalized  frequency  =  6/32 


Histognam  of  Cohenency  at  normalized  frequecy  =  V/32 
Mean=0 .  &6690  ,  Var=0 . 30906a-04 ,  assumed  model  order  =  1 


0.68  0.68  0.68  0.68  0.69  0.69  0.69  0.70  0.70  0.70 

Bin  Ceniers 

esiimation,  assumed  model  order  =  1,  (model  1),  data  length  IQl 
Histogram  of  estimated  SMC  at  normalized  frequency  =  7/32 


Histogram  of  Coherency  at  normalized  frequecy  =  8/32 
Mean=0. 54264  ,  Var=0. 26461e-04,  assumed  model  order  =  1 


0.53  0.53  0.53  0.54  0,54  0.54  0.55  0.55  0.55  0.56 

Bin  Cen+ers 

Fi£.__fl._r98  :  SMC  estimation,  assumed  model  order  =  1,  (model  1),  data  length  1G( 
Histogram  of  estimated  SMC  ot  normalized  frequency  =  8/32 


Hisiognam  of  Coherency  at  normalized  frequecy  =  *9/32 
Mean=0. 41 192  ,  Var=0 . 21366a-04 ,  assumsd  model  order  =  1 


Bin  Ceirlers 

El3.’__8.-_1.?9.  •  SMC  estimation,  assumed  model  order  =  1,  (model  1),  data  length  1024  pcir.i 
Histogram  of  estimated  SMC  ot  normalized  frequency  =  9/32 


Histogram  of  Coherency  at  normalized  frequecy  =  10/32 
MeQn=0 . 301 17 ,  Var=0 .  13661e-04 ,  ossumed  model  order  =  1 


Bln  Cen-ters 

‘  SMC  estimation,  assumed  model  order  =  1,  (model  1),  data  length  1024  point 
Histogram  of  estimated  SMC  at  normalized  frequency  =  10/32 


Cbservotions/Bin  Obsendiotis/Bin 


Hlsiognam  of  Cohor-ency  at  normalized  froquocy  =  11/32 
Mean=0 . 21074  ,  Vor=0 . 80343e-05,  assumed  model  order  =  1 


Bln  Ceniers 

Fig.  R.201  :  SMC  estimation,  assumed  model  order  =  1,  (model  1),  data  length  1024  point' 
Histogram  of  estimated  SMC  at  normolized  frequency  =  11/32 

Histogram  of  Coherency  at  normalized  frequecy  =  12/32 
Mean=0, 13935 ,  Var=0. 43223e- 05,  assumed  model  order  =  1 


Bln  Ceniers 

estimation,  assumed  model  order  =  1,  (model  1),  data  length  1024  point; 
Histogram  of  estimated  SMC  at  normalized  frequency  =  12/32 


Cbservntians/Bin  CbsffVDtiore/Bin 


His-togram  of  Coherency  at  normalized  frequecy  =  13/32 
Mean=0 . 65290e-01  .  Var=0 . 23633b-05  ,  assumed  model  order  =  1 


Bin  Cen-ters 

Fig.  P.203  I  SMC  estimation,  assumed  model  order  =  1,  (model  1),  data  length  1024  poirti 
Histogram  of  estimated  SMC  at  normalized  frequency  =  13/32 


Histogram  of  Coherency  at  normalized  frequecy  =  14/32 
Meon=0. 46992e-01 .  Var=0. 15520e-05,  assumed  model  order  =  1 


Bln  Centers 

Ei9.-_6.-_?®'!  '  SMC  estimation,  assumed  model  order  =  1,  (model  1),  data  length  1024  print: 
Histogram  of  estimated  SMC  at  normalized  frequency  =  14/32 


Hlsiognam  of  Coherency  ot  normolJzed  frequecy  =  15/32 
Mean=0. 23222e-01  ,  \'ar=0. 1328Be-05,  assumed  model  order  =  1 


Fig.  P.205  :  SMC  estimation,  assumed  model  order  =  1,  (model  1),  data  length  1024  pc 
Histogram  of  estimated  SMC  at  normalized  frequency  =  15/32 


Histogram  of  Coherency  at  normalized  frequecy  =  16/32 
Mean=0. 13179e-01 ,  Var=0. 13042e“05,  assumed  model  order  =  1 


Bin  Cen+ers 

‘  estimation,  assumed  model  order  =  1,  (model  1),  data  length  1024  pc 
Histogram  of  estimated  SMC  at  normalized  frequency  =  16/32 


Normalized  Frequency 


Fig.  PI.  207  :  SMC  esiimation,  assumed  model  order  =  1, 
data  length  256  points 


SMC  esiimation,  general  method,  data  length  256>  point 


0.0  0.1  0.2  0.3  0.4  0.5 

Normalized  Frequency 

Ei9.-__6.-_20Z- 1  ••  SMC  estimation,  assumed  model  order  =  1, 

data  length  256  points,  Average  estimoted  SMC  over  100  runs 


Normalized  sigmo**2  cf  estimated  S.M.C 


0.51  0.52  0.53  0.54  0.55  0.5e  0.57  0.5S  0.59  0.59 

Bin  Cenler-s 

:  SMC  estimation,  ossumed  model  order  =  1,  (model  1),  data  length  256  poin 
Histogram  of  estimated  SMC  at  normalized  frequency  =  2/32 


Cbserwticns/Bin  Cbsavitiore/Bin 


fi.118 


Histogram  of  Coherency  at  normalized  frequecy  =  3/32 
Mean=0 . 6'74 12  ,  Var=0 . 93503e- 04  ,  assumed  model  order  =  1 


Bln  Centers 

Fig.__R._2y^  :  SMC  estimation,  assumed  model  order  =  1,  (model  1),  data  length  256  poiris, 
Histogram  of  estimated  SMC  at  normalized  frequency  =  3/32 


Histogram  of  Coherency  at  normalized  frequecy  =  4/32 
Mean=0. 99501 ,  Var=0. 66641e-05.  assumed  model  order  =  1 


Bin  Centers 

Fi£.__R._21_2  :  SMC  estimation,  assumed  model  order  =  1,  (model  1),  data  length  256  poinis, 
Histogram  of  estimated  SMC  at  normalized  frequency  =  4/32 


Cbservotions/Bin  Cbserwtjons/Bin 


ft. 119 


Hlsiognam  of  Coherency  at  normalized  frequecy  =  B/32 
Mean=0.  95776 ,  Var=0. 29676s -04,  assumed  model  order  =  1 


Bin  Ceniers 

Fia-_£  ._2K3  :  SMC  estimation,  ossumed  model  order  =  1,  (model  1),  dato  length  255  points, 
Histogram  of  estimated  SMC  at  normalized  frequency  =  5/32 


Histogram  of  Coherency  at  normalized  frequecy  =  6/32 
Mean=0 . 63742 ,  Var=0 . 66927e- 04,  assumed  model  order  =  1 


Bln  Ceniers 

Fiq.  fi.214  ;  SMC  estimation,  ossumed  model  order  =  1,  (model  1),  doto  length  256  points, 
Histogram  of  estimated  SMC  at  normalized  frequency  =  6/32 


Histognam  of  Cohensncy  at  normalized  frequecy  =  7/32 
Mean=0 .  &6919  ,  Var=0 .  1 1725e-03,  assumed  model  order  =  1 


Bin  Cen-ters 

Fi£.__fi^._2^5  :  SMC  es-timation,  assumed  model  order  =  1,  (model  1),  data  length  256  points 
Histogram  of  estimated  SMC  at  normalized  frequency  =  7/32 


Histogram  of  Coherency  at  normalized  frequecy  =  8/32 
Mean=0. 54305 .  Var=0 .  10916a-03,  assumed  model  order  =  1 


Bln  Cen-ters 

Fi2^.__ft._21_6  !  SMC  estimation,  assumed  model  order  =  1,  (model  1),  data  length  256  points 
Histogram  of  estimated  SMC  at  normalized  frequency  =  8/32 


Cbservotions/Bin  Cbsavntiors/Bin 


His-togram  of  Cohenency  at  normalized  froquecy  =  9/32 
Mean=0 .  41 233  ,  Var=0 .  830‘99e-04 ,  (assumed  model  order  =  1 


Fig.  P.217  :  SMC  estimation,  assumed  model  order  =  1,  (model  1),  doto  length  256  poin:c, 
Histogram  of  estimated  SMC  at  normalized  frequency  =  902 


Histogram  of  Coherency  at  normalized  frequecy  =  10/32 
Mean=0 . 30153 ,  Var=0 . 55094e-04 ,  assumed  model  order  =  1 


Bln  Cen-ters 

Fi£.__fl._21t3  :  SMC  estimation,  assumed  model  order  =  1,  (model  1),  data  length  256  points, 
Histogram  of  estimated  SMC  at  normalized  frequency  =  10/32 


Observoticns/Bin  Cbsavitjons/Bin 


fi.l22 


Hisiogram  of  Coherency  of  normalized  frequecy  =  11/32 
Mean=0 . 21 102  ,  Var=0 . 32998s-04 ,  assumed  modal  order  =  1 


BJn  Cerrters 

Fi§_.__P._2]^9  :  SMC  estimation,  assumed  model  order  =  1,  (model  1),  doto  length  256  poii-'ts, 
Histogram  of  estimated  SMC  ot  normolized  frequency  =  11/32 


Hisfogram  of  Coherency  at  normalized  frequecy  =  12/32 
Mean=0. 13954 ,  Var=0,  16611e-04,  assumed  model  order  =  1 


Bln  Cen-ters 

Ei9.-__6.'_220  :  SMC  estimation,  assumed  model  order  =  1,  (model  1),  doto  length  256  points, 
Histogram  of  estimated  SMC  at  normalized  frequency  =  12/32 


Hlstognam  of  Coherency  of  normollzed  frequecy  =  13/32 
Mean=0. 65412e-01  ,  Var=0 .  10732e-04 ,  assumed  model  order  =  1 
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Bin  Ceniers 

Fig.  fl.221  !  SMC  esfimaiion,  assumed  model  order  =  1,  (model  1),  dofo  length  256  pcin: 
Histogram  of  estimated  SMC  at  normalized  frequency  =  13/32 


Hisfogrom  of  Coherency  of  normalized  frequecy  =  14/32 
Mean=0 . 47057e-0 1 ,  Var=0 . 72120e-05 ,  assumed  model  order  =  1 
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Bin  Cerrters 

Fig^.__fi._222  :  SMC  estimation,  assumed  model  order  =  1,  (model  1),  data  length  256  point 
Histogram  of  estimated  SMC  at  normalized  frequency  =  14/32 


CbssTvoticre/Sin  Cbservotions/Sin 


Histogram  of  Coherency  at  normalized  frequecy  =  15/32 
Mean=0 . 23250e-0 1 .  Var=0 . 60335e-05 ,  assumed  model  order  =  1 
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Bln  Cen-ters 

Fig.  P.223  !  SMC  estimation,  assumed  model  order  =  1,  (model  1),  data  length  256 
Histogram  of  estimated  SMC  at  normalized  frequency  =  15/32 


Histogram  of  Coherency  at  normalized  frequecy  =  16/32 
MeQn=0 .  13190e-0 1 .  Var=0 . 57363e-05 ,  assumed  model  order  =  1 
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Bln  Ceniers 

estimation,  assumed  model  order  =  1,  (model  1),  doto  length  255 
Histogram  of  estimated  SMC  at  normalized  frequency  =  16/32 


Normalized  Frequency 


Fig.  R.225  :  SMC  estimation,  assumed  model  order  =  1, 
data  length  128  points 


SMC  esiimatiion ,  general  method,  data  length  128  point 
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Fig^.__fl._225. 1  :  SMC  estimotion,  ossumed  model  order  =  1, 

data  length  128  points,  fiveroge  estimated  SMC  ouer  100  runs 
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Fig.  P.226  !  Squored  Siandaru 
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Fig.  fl. 226.1:  Normalized  Squi 

dala  lei 


Observoticns/Bin  Observotiore/Bin 


His+ognam  of  Coherency  of  normalized  frequecy  =  1/32 
Mean=0.  Var=0. 44390e-03,  assumed  model  order  =  1 


Bin  Cen-ters 

Fig.  R.227  :  SMC  estimation,  assumed  model  order  =  1,  (model  1),  data  length  128  poi". 
Histogram  of  estimated  SMC  at  normalized  frequency  =  1/32 


Hisfogrom  of  Coherency  of  normalized  frequecy  =  2/32 
Mean=0. 55994  ,  Var=0 . 33466e- 03,  assumed  model  order  =  1 


Bln  Ceniers 

Fig.  fi.228  :  SMC  estimation,  assumed  model  order  =  1,  (model  1),  data  length  128  pcin 
Histogram  of  estimated  SMC  at  normalized  frequency  =  2/32 


Histogram  of  Cohsrsncy  at  normal  Izsd  frsquscy  =  3/32 
Mean=0. 67378,  Var=0.  17577e-03.  assumed  model  order  =  1 
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esifmaiion,  assumed  model  order  =  1,  (model  1),  doio  length  128  points 
Histogrom  of  estimated  SMC  at  normalized  frequency  =  3/32 


Histogram  of  Coherency  at  normalized  frequecy  =  4/32 
Mean=0. 99367 ,  Var=0 .  13844e- 04,  assumed  model  order  =  1 
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Bin  Centers 

estimation,  assumed  model  order  =  1,  (model  1),  data  length  128  poin 
Histogram  of  estimated  SMC  at  normalized  frequency  =  4/32 


fl.l30 


Hlsiognam  of  Cohenency  at  normalized  frequecy  =  7/32 
Mean=0. 68725  .  Var=0. 2601 4e- 03,  assumed  model  order  =  1 


Bln  Cen-ters 


Fig.  P.233  :  SMC  es-timation,  assumed  model  order  =  1,  (model  1),  doio  length  125  points, 
Histogram  of  ei-timated  SMC  at  normalized  frequency  =  7/32 


Histogram  of  Coherency  at  normalized  frequecy  =  8/32 
Mean=0 . 54132 ,  Var=0 . 25424e-03,  assumed  model  order  =  1 


Bln  Centers 

Fig.  R.254  :  SMC  estimation,  assumed  model  order  =  1,  (model  1),  data  length  123  points, 
Histogram  of  estimated  SMC  at  normalized  frequency  =  8/32 


Histogram  oT  Cohensncy  at  normalized  frequecy  =  9/32 
Mean=0 . 410t>8  ,  VQr="0 .  16925a-0-3,  assumed  model  order  =  1 


Bin  Centers 

Fig.  R.235  s  SMC  estimation,  assumed  model  order  =  1,  (model  1),  data  length  128 
Histogram  of  estimated  SMC  at  normalized  frequency  =  9/32 

Histogram  of  Coherency  at  normalized  frequecy  =  10/32 
MeQn=0. 30036 ,  Var=0 .  12236e-03,  assumed  model  order  =  1 


Bln  Centers 

Fig.  R.236  :  SMC  estimation,  assumed  model  order  =  1,  (model  1),  data  length  126 
Histogram  of  estimated  SMC  at  normalized  frequency  =  10/32 


Hisiogr-am  of  Coherency  at  normalized  frequecy  =  11/32 
Mean=0 . 2101 1  ,  Van=0 . 709420 -04  ,  assumed  model  order  =  1 


Fig.  R.237  :  SMC  estimaiion,  assumed  model  order  =  1,  (model  1),  data  length  128  po 
Histogram  of  estimated  SMC  at  normalized  frequency  =  11/32 


Histogram  of  Coherency  at  normalized  frequecy  =  12/32 
Mean=0 .  13666 ,  Var=0 . 36325e- 04 ,  assumed  model  order  =  1 


Bln  Cen+ers 

estimation,  assumed  model  order  =  1,  (model  1),  data  length  128  po 
Histogram  of  estimated  SMC  at  normalized  frequency  =  12/32 


Cbser-iitiorB/Bin  Cbsavitjons/Bin 


fi.l33 


Histogram  of  Cohsrsncy  at  normal Izsd  frsquscy  =  13/32 
Mean=0. 64900e-01  ,  Var=0 . 21 1090-04  ,  assumed  model  order  =  1 


Bln  Centers 

:  SMC  esifmation,  assumed  model  order  =  1,  (model  1),  doio  length  128  poiris, 
Histogram  of  estimated  SMC  at  normalized  frequency  =  13/32 


Histogram  of  Coherency  at  normalized  frequecy  =  14/32 
Meon=0 . 46672e-0 1 ,  Var=0 .  13930e“04 ,  assumed  model  order  =  1 


Bln  Centers 

Fiq.  fl.240  !  SMC  estimation,  assumed  model  order  =  1,  (model  1),  doto  length  128  points, 
Histogram  of  estimated  SMC  at  normalized  frequency  =  14/32 


Observotions/Bin  Obsavjtiora/Bin 


His-fcogram  of  Coherency  of  normolized  frequecy  =  15/32 
Mean=0 . 22945s-01  .  Var=0 .  1 19020-04  ,  assumed  model  order  =  1 


Fig.  Pi. 241  :  SMC  esfima+ion,  assumod  Tcdel  order  =  1,  (model  1),  doto  length  12S 
Histogram  of  estimated  SMC  at  normalized  frequency  =  15/32 

Histogrom  of  Coherency  at:  normalized  frequecy  =  16/32 
Mean=0.  12919e-01 .  Var=0, 11651e-04,  assumed  model  order  =  1 


Bln  Cen-ters 

Fig.  fi.242  !  SMC  estimation,  assumed  model  order  =  1,  (model  1),  data  length  123 
Histogram  of  estimated  SMC  at  normalized  frequency  =  16/32 
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SMC  esilmaxlon ,  (model  1),  data  length  512  points 
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Fig.  R.244  :  Squared  Standard  Deviation,  of  estimated  SMC, 

assumed  model  order=  4,  100  runs 

SMC  estimation,  (model  1),  data  length  512  points 
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Fig.  R.244. 1:  Normalized  Squored  Standard  Deviation  of  estimoted  SMC, 

assumed  model  order=  4,  100  runs 


fi.l37 


Hisiognam  of  Cohorsncy  at  normalized  frequecy  —  1/32 
Model  1,  MBan=0 . 20405 ,  Var=0. 4&630e- 02,  assumed  model  order  =  4 


Fig.  R.245  s  SMC  estimation,  assumed  iTiOubI  order  =  4,  (model  1),  data  length  512  points. 
Histogram  of  estimated  SMC  at  normalized  frequency  =  1/32 


Histogram  of  Coherency  at  normalized  frequecy  —  2/32 
Model  1,  Mean=0. 56763,  Var=0, 356070-02,  assumed  model  order  =  4 


BJri  Cen'ters 

Fig.  F).246  :  SMC  estimation,  assumed  model  order  =  4,  (model  1),  data  length  512  coinip 
Histogram  of  estimated  SMC  ot  normalized  frequency  =  2/32 


Observotians/Bin  Oteovitjans/Bin 


fi.l38 


Hls+ognam  of  Coherency  oi:  normolized  frequecy  =  3/32 
Model  1,  Mean=0 •  BTVIB ,  Var=0.  490&7©-03 ,  assumed  model  order  =  4 


Bln  Ceniers 

Fig,  fl.247  z  SMC  esiimation,  assumed  model  order  =  4,  (model  1),  dcrto  length  512  points, 
Histogram  of  estimated  SMC  at  normalized  frequency  =  3/32 


His-togrom  of  Coherency  at  normalized  frequecy  =  4/32 
Model  1,  Mean=0. 99612,  Var=0.  16403e-05,  assumed  model  order  =  4 


B/ri  Ceirtors 


Fig.  R.248  ;  SMC  estimation,  assumed  model  order  =  4,  (model  1),  data  length  512  points. 
Histogram  of  estimated  SMC  at  normalized  frequency  =  4/32 


His'tognam  of  Coherency  of  normollzed  frequecy  “  7/32 
Model  1,  Mean=0. 66212,  Var=0. 2230'7e- 02,  oseumed  model  order  = 


.0 


Bln  Cen-tere 

Fig.  R.251  :  SMC  es-timaiion,  assumed  model  order  =  4,  (model  1),  doio  length  512  points, 
Histogram  of  estimated  SMC  at  normalized  frequency  =  7/32 


His-fcogrom  of  Coherency  at:  normalized  frequecy  =  6/32 
Model  1,  Mean=0. 53304 ,  Var=0.  44609e- 02,  assumed  model  order  =  4 


Bln  Centers 

Fi3^.__H._^2  !  SMC  estimation,  assumed  model  order  =  A,  (model  1),  data  length  512  points, 
Histogram  of  estimated  SMC  at  normalized  frequency  =  8/32 


Histogram  of  Cohsrsncy  at  normalized  frequecy  =  9/32 
Model  1,  Mean=0  •  405&7 ,  Var=0 .  50260e-02 ,  assumed  model  order  = 
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Fig,  fl.253  :  SMC  estimation,  assumed  model  order  “  4,  (model  1),  data  length  512  pcirits. 
Histogram  of  estimated  SMC  at  normalized  frequency  =  9/32 


Histogram  af  Coherency  at  normalized  frequecy  =  10/32 
Model  1,  Mean=0. 30212,  Var=0. 40191b-02,  ossumed  model  order  =  4 
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Bln  Cen-ters 

Fig,  fl.254  !  SMC  estimation,  assumed  model  order  =  4,  (model  1),  data  length  512  points. 
Histogram  of  estimated  SMC  at  normalized  frequency  =  10/32 


12 


Histogram  of  Cohsrsncy  at  normal Izsd  frsquscy  =  11/32 
Model  1,  Mean=0. 21657,  Var=0. 34507b-02,  assumed  model  order  =  4 


Bln  Cen"ler3 

Fig.  R.255  i  SMC  estimation,  assumed  model  order  =  4,  (model  1),  doto  length  512  points, 
Histogram  of  estimated  SMC  at  normalized  frequency  =  11/32 


Histogram  of  Coherency  at  normalized  frequecy  =  12/32 
Model  1,  Mean=0. 15253,  Var=0. 35224e- 02,  ossumed  model  order  =  4 


Bln  Cen-tors 

Fig.  R.256  i  SMC  estimation,  assumed  model  order  =  4,  (model  1),  data  length  512  points, 
Histogram  of  estimated  SMC  at  normalized  frequency  =  12/32 


HJtsiogram  of  Coherency  at  normalized  frequecy  =  13/32 
Model  1,  Maan=0. 10217 ,  Var=0. 26976e-02,  assumed  modal  order  = 
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Fig.  R.257  t  SMC  es-timation,  assumed  model  order  =  4,  (model  1),  data  length  512  points 
Histogram  of  estimated  SMC  at  normalized  frequency  =  13/32 


Histogram  of  Coherency  at  normalized  frequecy  =  14/32 

Model  1,  MBan=0. 64666e-01 ,  Var=0.  16151b-02,  assumed  model  order  =  4 
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Bln  Ceniers 

:  SMC  estimation,  assumed  model  order  =  4,  (model  1),  data  length  512  points 
Histogram  of  estimated  SMC  at  normalized  frequency  =  14/32 


(bservotiofiyBin  Cbeendions/Bin 


fi.l44 


Histogram  of  Cohorency  at  normal fzsd  froquecy  =  15/32 

Model  1,  MBan=0. 39360e- 01 ,  Var=0. 92695e- 03,  assumed  model  order  =  4 


Bin  Ceniers 

Fig.  P.259  :  SMC  es-timation,  aBsumed  model  or-der  =  4,  (model  1),  doio  length  512  points, 
Histogram  of  estimated  SMC  at  normalized  frequency  =  15/32 


Histogram  of  Cohsrency  at  normalized  frequecy  =  16/32 

Model  1,  MBan=0. 274Q4e-01 ,  Var=0.  10265e-02,  assumed  model  order  =  4 


Bln  Cen-lers 

Fig.  R.260  t  SMC  estimation,  assumed  model  order  =  4,  (model  1),  data  length  512  points, 
Histogram  of  estimated  SMC  at  normalized  frequency  =  16/32 


